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ground and large-scale structure on the models of primordial non-Gaussianity. We study 
the generalized local model of non-Gaussianity, where the parameter /nl is promoted to 
a function of scale, and present the principal component analysis applicable to an arbi- 
trary form of /nl(&)- We emphasize the complementarity between the CMB and LSS by 
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1. Introduction 

There has recently been a surge in interest to study departures in the distribution of primor- 
dial density fluctuations from the random Gaussian case predicted by standard inflationary 
models. The reason for this renewed interest lies in the fact that any observable departures 
from Gaussianity would essentially rule out the standard single-field, slow-roll inflationary 
picture, pointing to a more complicated dynamics during the epoch of inflation (see e.g. 
[1, 2] for reviews). 
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It is therefore important to consider how one could parametrize primordial non-Gaussianity. 
A much-studied model of primordial non-Gaussianity is the local (or squeezed) model, 
which characterizes non-Gaussianity through a single parameter /nl [3, 4, 5] 

$(X) = 4> G {X) + fNL(Mx? ~ {4>g{x?)). (1.1) 

Here, $ denotes the primordial curvature perturbations (Bardeen's gauge-invariant poten- 
tial), <I)g( x ) is a Gaussian random field, and the constant /nl is the parameter describing 
deviations from Gaussianity. The local model has been much studied for its simplicity - 
it contains the first two terms of the most general local form of non-Gaussianity [6]. In a 
recent paper ([7]; hereafter BHK11), we introduced a generalization of this model to one 
where, in Fourier space, /nl = /nl(&) is a function of scale 

$(fc) = <j> G {k) + f NL (k) J —<j> G (k/)<l> G (k - k'). (1.2) 

This is a natural extension of the popular 'local' model 1 /nl = const, and it is physically 
well-motivated; it can describe inflationary scenarios with multiple light fields, one of which 
is responsible for the generation of curvature perturbations [9, 10, 11]. The bispectrum in 
this generalized local model is 

B4h,k 2 , k 3 ) = 2[f NL (ki)P (f) (k 2 )P (l> (k 3 ) + perm.], (1.3) 

where Pa is the power spectrum of potential fluctuations. This reduces to the familiar 
expression B(k±, k 2 , ks) = 2fj<n J (P ( f ) (ki)P ( f l (k 2 ) +perm.) when /nl is a constant. 

To parametrize this model while retaining its full generality, it is convenient to consider 
a parametrization in piecewise-constant bins in wavenumber: 

/NL = fNL(h) (1.4) 

where each /^j L is the value of /nl(^) in the i-th wavenumber bin. In BHK11, we used this 
parametrization to project errors on /nl(&) from a hypothetical Stage III galaxy survey. 
As in BHK11, we adopt 20 bins in wavenumber distributed uniformly in log(fc), which is 
easily sufficient to obtain the best-measured principal components accurately. 

In this work, we perform an analysis that is similar in spirit to that in BHK11, but 
extended in several respects. First of all, we develop formalism and work out forecasts for 
how well the CMB, in particular Planck [12], can measure /nl(^)- We combine this with 
the LSS forecasts updated to reflect three specific galaxy surveys. Having done that, we 
obtain a clearer picture of where we can expect good constraints on non-Gaussianity in 
/c-space. Finally, we project our forecasts to the specific, power-law model in wavenumber, 
and thus clarify at which wavenumber the CMB, LSS and the combined surveys best 
determine non-Gaussianity. Therefore, this work complements not only BHK11 and studies 

1 Other models also display scale dependence of primordial non-Gaussianity; for example, the Dirac-Born- 
Infeld braneworld theory typically leads to scale-dependent equilateral /nl, which has been constrained in 
Ref. [8]. 
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that forecasted errors for (occasionally slightly different) models of scale-dependent non- 
Gaussianity [13, 10, 14], but also many previous forecasts for future constraints on constant 
/ NL [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 14, 27, 28]. 

The structure of this paper is as follows: in Sec. 2, we briefly review the main result 
from BHK11 - the signature of the generalized local model on large-scale structure through 
halo bias - and explore the effects of an additional term in the modeling of non-Gaussian 
bias first pointed out by Desjacques et al. [29]. In Sec. 3, we find the signature of the 
generalized local model on the CMB bispectrum, with particular emphasis on Planck. 
Finally, in Sec. 4 we combine the results for a set of joint constraints; we also perform a 
principal-component analysis, and project constraints on a power-law model of /nl(^)- We 
conclude in Sec. 5. Details of the computational work can be found in the Appendices. 



2. Generalized local model: signatures in large-scale structure 

In this section we first briefly review how the general local model of non-Gaussianity affects 
the bias of dark matter halos. We then present details of the LSS surveys that we will 
consider. 



2.1 Effect on the bias 

For the local non-Gaussian model from Eq. (1.1), the dark matter halo bias acquires scale 
dependence [30]: 

on tj2 

«*) = fa 4- A„ W = fa + / NL (fa - 1)4 a8(a)r " (t) V (2.1) 

where &o is the usual Gaussian bias (on large scales, where it is constant), 5 C ~ 1.686 is the 
collapse threshold, a is the scale factor, Q m is the matter density relative to the critical 
density, Hq is the Hubble constant, k is the wavenumber, T(k) is the transfer function, 
and g(a) is the growth suppression factor 2 . See also [31, 32, 33, 34, 18, 35, 36, 37, 38, 39, 
24, 40, 41, 42, 43, 44] who explored the effects of primordial non-Gaussianity on the bias 
of dark matter halos in great detail. 

In BHK11, we worked out the signature of the generalized local model, using the MLB 
formalism [45, 46, 31], and made forecasts for future galaxy surveys. It is not our intention 
to fully repeat all of our analysis from that paper, and we just quote essential results. The 
change in the bias Ab is given, in the generalized local model, by Eq. (3.16) from BHK11 



x / d^M R {k 2 ) 



hUk)^l + 2h L (k 2 ) 
p 4>(k) 



(2.2) 



2 The usual linear growth D(a), normalized to be equal to a in the matter-dominated epoch, is related 
to the suppression factor g(a) via D(a) — ag(a) / g(l), where g(a) is normalized to be equal to unity deep 
in the matter-dominated epoch. 
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where Mn(k) = [k 2 T(k)Wn(k)]/(HQQ m ) and Wji(k) is the Fourier transform of the top- 
hat filter with radius R. The derivatives with respect to piecewise constant parameters 
/Jj L are straightforward and given by Eq. (3.19) of our previous paper. 

However Eq. (2.2) is not entirely correct in describing the scale-dependent bias from 
a general NG model. Desjacques et al. ([29]; see also [41]) pointed out that the expression 
(2.1) is only correct in the high-peak, small-A; limit. An additional term is required for the 
exact expression: 



Ab(k) 



2F(k) 



M R {k)D{z) 



(b - 1)S C + 



rflnF(fc) 
din (Jr 



where 



F(k) 



8vr 2 



a 



R 



dhkfMRik^P^h) / d f iM R (k 2 ) 



/nl(*0 



+ 2h L {k 2 ) 



(2.3) 



.(2.4) 



The new term (second term in square parentheses in Eq. (2.3)) vanishes when the 
fiducial model for non-Gaussianity is /nl(&) = const, but it becomes relevant for truly 
scale-dependent models, including the piecewise-constant parametrization of /nl(&) from 
equation (1.4). Because we are expanding (taking Fisher derivatives) around the constant 
value of /nl(^) (30 or zero), we find very small though nonzero effect of this new term 
describing halo bias. See Appendix C for details. 

2.2 Fisher matrix analysis: assumptions and survey specifications 

We are interested in making forecasts for constraints on non-Gaussianity from future galaxy 
surveys, and for this we employ the standard Fisher matrix formalism. 

For measurements of the power spectrum of dark matter halos, the Fisher matrix F is 

[47] 

/■*»- dP h (k,z m )dP h (k,z m ) 1 k 2 dk 

femin 



hLSS 



dpi 



d Pj 



Ph{k,z m ) + 



1 



n 



(27T) 



(2.5) 



where V m is the comoving volume of the m-th redshift bin, each redshift bin is centered 
on z m , and we have summed over all redshift bins. We adopt 3 k m \ n = 10 -4 h^ 1 Mpc, 
and we choose fc max as a function of z so that <7(7r/(2fc max ), z) = 0.5 [48], which leads to 
&max(-z = 0) « 0.1/iMpc -1 . Ph is the dark matter halo power spectrum, related to the true 
dark matter power spectrum P through 



P h (k) =b(k) 2 P(k), 



(2.6) 



where each quantity implicitly also depends on redshift. Finally, pi are the parameters 
of interest; in our case, these are the f^, cosmological parameters, and the bias-related 



3 We found that the constraints are insensitive to the precise value of fc m i n for the fiducial /nl(A;) = 30 
(or any sufficiently nonzero value), since the most constraining scales are intermediate between fc m i n and 



fen 



reflecting the competition between larger noise and larger signal as one goes to lower k. However, for 



the fiducial value /nl(&) = the constraints indeed come from the largest available scales, and in that case 

We adopt fc m i n = VsurVey- 
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nuisance parameters listed below in (2.7). The minimal error in the i-th cosmological 
parameter is, by the Cramer-Rao inequality, o~(pi) ^ \J (F^ 1 )a. 

All of the results cited in this section, as well as LSS survey projections elsewhere in this 
paper, assume the following survey properties, modeled on BigBOSS [49], unless explicitly 
stated otherwise. We adopt the fiducial value /nl = 30, chosen to roughly correspond to 
the maximum-likelihood value favored by current CMB data [50]. The fiducial cosmological 
model is the standard ACDM model with Hubble's constant Ho; physical dark matter and 
baryon densities J7 c d m ^ 2 and Q^h 2 ; equation of state of dark energy w; the log of the scalar 
amplitude of the matter power spectrum, log A s ; and the spectral index of the matter power 
spectrum, n s . Fiducial values of these parameters correspond to their best-fit WMAP7 
values [50]. We also added the forecasted cosmological parameter constraints from the CMB 
experiment Planck by adding its Fisher matrix as a prior (W. Hu, private communication). 
Note that the CMB prior does not include CMB constraints on non-Gaussianity, so that we 
are not double-counting the latter constraints. Next, we include twenty piecewise-constant 
non-Gaussianity parameters /nl(&i) = /nl with i = 0,1,..., 19. Finally, we include a 
Gaussian bias parameter in our Fisher matrix, bo(z), for each of the 44 redshift bins over 
the range 0.1 < z < 4.5 (which is discussed further below). The full list of parameters that 
we have in the BigBOSS Fisher matrix is 

{ Pi } = {H , n cdm h 2 , Q h h 2 , w, log A s , n s , bf, / NL , . . . , / N 9 L } (2.7) 

BigBOSS will utilize three different tracers: LRGs, ELGs, and QSOs. The fiducial values 
for bo(z = 0) are different for each of the three tracers, as are the number densities. To 
account for this, we calculated the Fisher matrix for each tracer, then added all three 
together for the final BigBOSS Fisher matrix. The fiducial values for bo(z = 0) were 1.7 
for the LRGs, 0.84 for the ELGs, and 1.2 for the QSOs [51]. In each case, we assume the 
simple scaling of bo with redshift, bo(z) = bo(z = 0)/D(z). 

In order to simplify the analysis, and in light of the uncertainties of the distribution 
of observed LSS tracers' masses, we assume a fixed halo mass of I^Mq/Ji. Since we 
marginalize over bo in each redshift bin, we effectively delete much information about the 
mass of the tracers. Essentially, we utilize information about the redshift- and wavenumber- 
dependence of bias, but avoid — at least for now — using information about the masses, 
since accurate masses of LSS tracers are typically very difficult to obtain, except for galaxy 
clusters. 

We assume that BigBOSS will cover 14,000 square degrees. The redshift ranges for 
the LRGs, ELGs and QSOs are zg%* = 1.2, z^G = L8 and z QSO = 4 _ g . thg forecasted 

number density in each redshift bin of course widely vary, and in particular the QSOs have 
much lower typical number densities per redshift bin than the other tracers [51]. We split 
the survey into the aforementioned 44 redshift bins out to z max = 4.5; the total volume of 
the survey is therefore V tot = 230 (/i^Gpc) 3 . The errors in the cosmological parameters 
vary, in the cosmic variance limit, as V to ^ 2 . 

In addition to presenting our results with the fiducial BigBOSS survey (14,000 sq. deg., 
-Zmax = 4.5), we also forecast constraints for the Dark Energy Survey (DES [52]), which 
will cover 5,000 sq. deg. with 2; max = 1-0, split into 5 redshift bins. For simplicity's sake, 
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we assume that this survey will only see one tracer, with a fiducial Gaussian bias bo(z = 
0) = 2.0 and a number density n = 2 x 10~ 4 (/iMpc" 1 ) 3 (independent of redshift). Note 
that the actual redshift distribution and number density of the DES tracers will depend 
on the spectroscopic followup to the main survey and the objects that it targets, details of 
which are somewhat uncertain at this time and hence our approximate assumptions. The 
rest of the Fisher matrix formalism for this survey (cosmological parameters, etc.) is the 
same as what we used for the BigBOSS Fisher matrix. 

2.3 The effect of the fiducial value on constraints 

Looking back at Equation (2.5), we see that the fiducial /nl enters through the bias, by 
way of Pfi- Assuming Ph(k) S> 1/n (a reasonable assumption at large angular scales where 
non-Gaussianity constraints largely come from and where shot noise is negligible), we find 
that the Fisher matrix element corresponding to /nl = const is 

Thus, the expression on the right-hand side will, in general, be dependent on the choice of 
fiducial /nl- Since |A6(/c)| blows up at small k, in that regime we have: 

Ab(k) \ 2 _ J_ 
'nl 



/nl (bo + Ab(k)) J A 2 (2 ' 9) 



At large k, Ab(k) goes to 0, taking the entire expression with it. Thus, the integral is 
dominated by the contribution at low k, meaning we should expect a maximal Fisher 
matrix element around a fiducial /nl = 0. And indeed, that is what we see in Figure 1: 
the projected constraints on /nl from a given sky survey depend strongly on the fiducial 
value chosen, with the tightest constraints at /nl = 0. 

We note that, as shown in the following section, the Fisher matrix is independent of 
the fiducial /nl value for the CMB constraints for our piecewise-constant parameterization. 



3. Generalized local model: signatures in the CMB 

Traditionally, the best constraints on non-Gaussianity have come from the CMB. This is 
done almost exclusively through estimators involving the N-point correlation functions for 
N > 2 and their Fourier transforms, the polyspectra. Most emphasis has been on the 
N = 3 case, or the bispectrum of temperature fluctuations in the CMB, if only because 
of its relative computational simplicity. The well-known general expression for the CMB 
bispectrum, re-derived in Appendix A, is 

"SS» - (§) 3 y MHEE Q £ *) J ^ «*, ,1*3 

/>oo 

x B^(ki,k 2 ,k3)e ei (ki)tUk2)f i3 (h) / r 2 dr j h (kir)u 2 (k 2 r)ji 3 (k 3 r) (3.1) 

Jo 
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Figure 1: A more detailed look at how the choice of fiducial /nl affects the projected constraints 
on constant /nl from a future galaxy survey. See text for analytic explanation for why results are 
the best at fiducial value of /nl = 0. 



where the expression in angular parentheses is the Wigner-3j symbol, B§ is the curvature 
bispectrum, and t( are the radiation transfer functions. 

In principle, we can use this to find the Fisher matrix F^j for the CMB bispectrum and 
thus forecast how well the CMB bispectrum can determine the cosmological parameters: 
[53, 54, 5, 55] 



-. flp>lmn BB Pqr 

F CMB _ r \^ \^ 1 un 'hMs (ri-l \ OD lxt-ih (o \ 

b H -^ky 2^ A dp . ^ih£ 3 )lmn, Pq r g ■ 13.2 ) 



Here, C is the covariance of the bispectra and pij are the parameters of interest. Af 1 ^ 2 ^ 3 
is a combinatoric term - equal to 6 when £\ = £2 = £3, 1 when l\ ^ £2 7^ £3, and 
2 otherwise [55]. The indices k and p,q,r run independently over all eight possible 
ordered triplets of temperature and polarization fields (TTT, TTE. . . EEE). The details of 
calculating and its derivatives are in Appendix A, while the details of calculating 

the bispectrum covariance C are in Appendix B. 

Equation (3.1) is a totally general result for the bispectrum of the CMB in terms of 



-7- 



the Bardeen curvature bispectrum; we have not picked any model of non-Gaussianity. But 
(3.1) is not useful without picking a form for B$(ki, ks). For the constant /nl case, we 
have the following Bardeen curvature bispectrum: 

B*(h, k 2 , k 3 ) = 2A2/ NL ^^r^i^p^ziy + Perm.^J (3.3) 

where A^ is the amplitude of the curvature power spectrum. Using Eqs. (3.2), (B.4), and 
(A. 24), we have the following expression for the CMB bispectrum Fisher information in 
the constant /nl case: 

1 (24 + 1) (2£ 2 + l) (24 + 1) (hhhV 1 



pCmb _ < a4 v v l (24 + i)(24 + i)(24 + i) /4 4 4\ 
/nl - Z_> a 4?r I o / a 



lmn,pqr 2<i^<i 2 <iz 



oo 



r 2 dr (af tl (r)P%(T)ft a (r)+j>erm. 



(3.4) 



oo 



dr (a^(r)/^ 2 (r)/3[ 3 (r) + perm. 



where ag and j3g are defined in equations (A. 21) and (A. 22). 

For the scale- dependent /nl(&) case from our generalized ansatz, things are somewhat 
more complicated. The Bardeen curvature bispectrum is: 

B^MM) = 2A^ ( fc 3-J-i L )^-K-i) +perm -) • (3 - 5) 

Using the piecewise-constant parametrization of /nl(^)j Eqs. (3.2), (B.4), and (A. 25) yield 
the following expression for the Fisher matrix of all the in the scale-dependent case, 
similar to Eq. (3.4): 



FT- 



CMB 



E E 

lmn,pqr 2<l 1 <l 2 <h 



l (24 + l)(24 + l)(24 + l)/4 



A, 



4tt 







A 



(C £i )ip(C £2 )jq{C £3 )kr 



r 2 dr (o?il{r)PT 2 (r)Pl{r) + perm. 



(3.6) 



r dr [ct^{r)f3 i2 (r)f3l 3 (r) + perm. 



Despite appearances, this is a relatively straightforward calculation to perform, and it takes 
roughly an hour (on four processors) for twenty parameters with ^ ma x ~ 

2000. 

We did not include other cosmological parameters in the CMB bispectrum Fisher ma- 
trix, since the bispectrum does not constrain them terribly well, while on the other hand the 
CMB power spectrum places very good constraints on the other cosmological parameters. 
In other words, non-Gaussianity estimates obtained using the CMB bispectrum would not 
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fc(h/Mpc) fe(h/Mpc) 
(a) BigBOSS (b) Planck 




k (h/Mpc) 

(c) Combined 

Figure 2: Constraints on the piecewise constant parameters in the generalized local model with 
the LSS (top left), CMB (top right), and LSS+CMB (bottom). All constraints are unmarginalized, 
in order to more clearly show the wavenumber-dependent sensitivity of the probes to primordial 
non-Gaussianity. The LSS constraints come from the power spectrum of galaxies, while the CMB 
constraints come from the bispectrum of temperature fluctuations. See text for details. For refer- 
ence, the green line is the constraint found for a constant /nl using the same assumptions. There 
are bins "missing" on the rightmost end of the Planck plot; those bins correspond to k- values too 
large to be probed when £ max = 2000, as it is here. 



be significantly affected by marginalizing over the other cosmological parameters within 
their allowed ranges, as explicitly shown by Ref. [56]. Therefore, it is a fair (and cer- 
tainly very helpful) approximation to think the CMB power spectrum and the bispectrum 
complementing each other by constraining the standard cosmological parameters, and the 
non-Gaussian parameters, respectively and separately. This has indeed been the approach 
in the literature (e.g. [5, 54]). 
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(c) Combined PCs 



Figure 3: The first four forecasted principal components of /nl(&) from BigBOSS, Planck, and 
BigBOSS+Planck, assuming the fiducial model fnh(k) = 30. The PCs eigenvectors e^\k) are 
ordered from the best-measured one (j = 0) to the worst-measured one (j = 19; not shown here) 
for the assumed fiducial survey. 



4. Results and Joint Constraints 
4.1 Forecasted constraints on the f^ L 

Figure 2 shows the (unmarginalized) constraints on the piecewise constant parameters 
/Jf L in the generalized local model from BigBOSS and Planck individually, as well as 
combined. Note that the two types of surveys have comparable constraints at the pivot 
wavenumber, but the pivot is at a larger wavenumber for BigBOSS. Away from the pivot, 
the Planck constraints are expected to be better than those from BigBOSS, but both 
rapidly deteriorate away from their respective pivots. Finally, the combined constraints are 
significantly helped by the lever arm in wavenumber when the two probes are combined, 
and this leads to better constraints across a wider range of scales. We will make these 
statements more quantitative below when we study the specific case where /nl(^) is a 
pure power law in k. 

The horizontal green curves in all panels in Fig. 2 show the accuracy in the constant 
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4 6 8 10 12 

Principal Component 



(b) Fiducial /nlW = 0, BigBOSS and Planck 



Figure 4: RMS error on each principal component for BigBOSS, Planck, and the two combined. 
Note that the BigBOSS errors are slightly smaller than those from Planck; in all cases combining the 
CMB and LSS decreases errors. Note too that for fiducial /nl(&) = (right panel) BigBOSS only 
constraints /nl(^) at the pivot point well (see Fig. 7), and hence the error in the best-determined 
principal component is noticeably smaller than errors in the other PCs. 



/nl, projected down from the principal components as described in BHK11. The 
accuracy achieved in /nl is 4.4 for Planck, 2.6 for BigBOSS, and 2.2 for the combined 
case. Recall also that our Fisher matrices for Planck - but not for BigBOSS - assume all 
cosmological parameters other than the f^ L are fixed (known). 

4.2 Principal Component Analysis 

Following BHK11, we now represent a general function /nl(&) i n terms of principal com- 
ponents (PCs). In this approach, the data determine which particular modes of /nl(&) are 
best or worst measured. The PCs also constitute a useful form of data compression, so that 
one can keep only a few of the best-measured modes to make inferences about the function 
/nl(&)- The PCs are weights in wavenumber with amplitudes that are uncorrelated by 
construction, and they are ordered from the best-measured (i = 0) to the worst-measured 
(i = 19) for the assumed fiducial survey. We follow the construction of the PCs following 
the formalism outlined in Appendix B of BHK11. We assume a total of 20 principal com- 
ponents distributed uniformly in log(10 -4 h Mpc" 1 ) < log(fe) < log(l /iMpc" 1 ), which is 
easily sufficient to make model-independent statements about /nl(&). 

Figure 3 shows the forecasted PCs of LSS and Planck separately and combined. Heuris- 
tically, the lowest principal component (PC0) serves to see how well we can find the devia- 
tion of /nl(&) & t its pivot (i.e. best-determined wavenumber) from the fiducial value. The 
higher PCs (PCI, PC2, etc) serve to probe the k-dependence of /nl- 

Figure 4 shows the l-cr errors on the PCs for BigBOSS, Planck, and the two combined. 
Note that the BigBOSS errors are slightly smaller than those from Planck (the DES er- 
rors, not shown here, are bigger than Planck's). Combining BigBOSS and Planck sharply 
decreases the errors. Note too that for fiducial /nl(&) — 0, BigBOSS only constraints 
/nl(&) a t the pivot point well (as shown below in Fig. 7), and hence the error in the 
best-determined principal component is noticeably smaller than errors in the other PCs. 

The relative strength of the LSS and CMB constraints at their respective pivot points 
strongly depend on two factors: volume of the LSS survey and, to a slightly lesser extent, 
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fiducial (i.e. true) value of /nl(&) (the CMB is not as sensitive to the fiducial value of 
/nl(&))- For example, for /nl(&) = 30 and DES we find that Planck constraints at the 
CMB pivot are stronger, while assuming BigBOSS survey we find that the LSS is slightly 
stronger at its own pivot. In addition, the CMB typically constrains /nl(&) ° ver a wider 
range of scales than the LSS. 

4.3 Projecting constraints on the power- law model of /nl(&) 

Once the Fisher matrix F has been obtained for the set of parameters /ml, it is quite simple 
to find the best possible constraints on the that could be obtained from a future galaxy 
redshift survey. By projecting this Fisher matrix into another basis, it is also possible to 
find the constraints on any arbitrary /nl(&) without calculating a new Fisher matrix from 
scratch. 

Here we will study the popular simple form of non-Gaussianity analogous to the con- 
ventional parameterization of the power spectrum [57, 58, 59, 60, 61, 13] 

/NL(fe) = fa f ¥ ) , (4-1) 

where fc* is an arbitrary fixed parameter, leaving /^ L and n/ NL as the parameters of interest 
in this model. The partial derivatives of our basis of with respect to these parameters 
are: 



Bf ^ ~ /nl f^V" 1 log I , (4.3) 



where ki is the k at the center of the zth /c-bin. Starting in a basis of 20 f^ L evenly spaced 
in log k, we project down to a basis of and nj NL in order to forecast constraints on 
the two new parameters. [Note that k* is an arbitrarily chosen parameter which differs in 
general from the true pivot k v - n where constraints on /nl(&) are the best. The choice of 
k* affects neither the constraints on /nl(&) nor the value of k p i v .] 

We can use the constraints on /^ L and nf NL to find constraints on /nl(&) as a whole, 
through the usual methods of error propagation: 

where C/^ L ,n /NL is the covariance between and n/ NL , and cr(f^ L ) 2 and cr(n/- NL ) 2 are 
their respective variances. Using this relation, and given some fiducial model of /nl(&)j 
we can plot the forecasted constraints on /nl(&) as a function of k. This is what we have 
done in Figure 5 for the Planck bispectrum, DES power spectrum, and the two combined 
(along with priors on cosmological parameters from the Planck power spectrum). 

Figure 6 is analogous to Figure 5, but shows BigBOSS and Planck constraints (rather 
than DES and Planck) for the fiducial value of /nl(&) = 30. Note that the forecasted 
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10* 10" 3 10" 2 10" 1 10° 



k (h/Mpc) 

Figure 5: Forecasted constraints on /nl(^) from several different data sets, assuming the power- 
law model of scale-dependent non-Gaussianity: /nl(^) = /NL(' c Apiv) n/NL > projecting down from 
the piecewise-constant basis. The red dashed line is the maximum k for which information was 
kept in the LSS Fisher matrix at z = 0. The LSS survey used for this forecast is based on DES. 



BigBOSS constraints are, very roughly, comparable to those from Planck (see also Table 
1), but are also very complementary to Planck since their best constraints are at a higher 
k. Our forecasted constraints on the accuracy of measuring the running with BigBOSS are 
in good agreement with forecasts for the Euclid survey in Ref. [14]. 

We also introduce the Figure of Merit (FoM( NG )) of non-Gaussianity. We defined it 
analogously to the Figure of Merit for dark energy ( [62] ; see also [63] ) as 

FoM( NG ) = (detF 2x2 ) 1 /2« 6 i 77r (4.5) 

^95 

where i*2x2 is the Fisher matrix projected down to the 2x2 space of and nf NL , 
and A95 is the area of the 95.4% confidence level ellipse in this space. Constraints on 
the FoM( NG ) are presented in Table 1, and show that combining of BigBOSS and Planck 
improves constraints by a factor of between two and five relative to these experiments 
alone. What is particularly encouraging is that future constraints will improve the recently 
obtained current constraints on the running of non-Gaussianity [64] by more than an order 
of magnitude. 
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lOOf 




k (h/Mpc) 



Figure 6: The same as Figure 5, but with survey parameters for large-scale structure based on 
BigBOSS. 



The constraints on /nl(^) from a large-scale structure survey are quite sensitive to 
the survey parameters. Unlike the constraints on /nl(&) from the CMB bispectrum, the 
forecasted constraints from LSS are also sensitive to the choice made for the fiducial model 



Projected errors c(/nl) an d <7( n f NL ), ana ^ * ne corresponding pivots 



Variable 


BigBOSS 


BigBOSS+Planck C e s 


Planck bispec 


BigBOSS+all Planck 


'(fa.) 


3.0 


2.6 


4.4 


2.2 


^/nl) 


0.12 


0.11 


0.29 


0.078 


FoM ( NG ) 


2.7 


3.4 


0.78 


5.8 


kpiv 


0.33 


0.35 


0.080 


0.24 



Table 1: Forecasted constraints on and n/ NL from BigBOSS, Planck, and combined data sets 
for two fiducial values of /nl(&)- Each column's numbers are for the pivot in that column; thus the 
errors in the two parameters are uncorrelated in each column. See text for survey specifications. 
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Projected errors (<7/^ L > & nf ) for different fiducial /nl(&) 





DES 


BigBOSS 


Planck 


Fiducial f NL (k) = 30 


(13, 1.0) 


(2.6, 0.11) 


(4.4, 0.29) 


Fiducial /nl(&) = 


(13, oo) 


(2.5, oo) 


(4.4, oo) 



Table 2: Forecasted constraints Cf* from different LSS surveys, assuming different fiducial models. 
Forecasted constraints from Planck are also shown for comparison. (All values of fi/ NL are equally 
likely in the second fiducial model, where f^ L = 0. ) 

of /nl(&) 5 as shown in Section 2.3. Forecasted constraints on f^ L and n/ NL for the DES and 
BigBOSS surveys, with two different fiducial models, are compared to forecasted constraints 
from Planck in Table 2 (note that all values of ^/ NL are equally likely for the fiducial 
model where f^ L = 0, and hence an infinite error on rif NL ). The scale at which the 
LSS gives the best constraint (the 'sweet spot') turns out to be slightly smaller than the 
maximum wavenumber assumed to be used by the survey, /c max - This is because the 
halo-bias integration over all the possible momentum space configurations in Eq. (2.2) has 
dominant contributions from small scales 4 , as we noted previously in [7]. Figure 7 shows the 
same case as Fig. 6, except for the fiducial value of /nl(^) = 0. Because nj NL is arbitrary 
for this fiducial value, constraints on /nl(^) are only good at a single, pivot wavenumber; 
this can also be seen by inspection of Eqs. (4.3) and (4.4). Even in this case (which, note, 
has measure zero in parameter space), we see that combination of BigBOSS and Planck is 
extremely beneficial 5 . 

5. Conclusions 

This paper focused on the ability of upcoming LSS and CMB surveys to probe more general 
models of primordial non-Gaussianity. We concentrated in particular on the generalized 
local model where the parameter /nl is promoted to an arbitrary function of scale /nl(^)- 
Our starting point were the piecewise constant parameters in k, constraints on which are 
shown in Fig. 2, and their principal components which are shown in Fig. 3 and constrained 
in Fig. 4. 

Comparison with theory is easiest, however, by using a simpler parametrization in 
terms of "running" of the spectral index, rif NL = din (k)/d In k. Using the two- 
parameter description of non-Gaussianity in terms of amplitude f^ L and running ny NL , 

4 We performed our bias calculations in the Lagrangian picture where the primordial fluctuations are 
linearly extrapolated to z = as usually done in the literature. For an alternative approach including the 
higher order corrections in the framework of the integrated perturbation theory, see Ref. [65] . 

5 While it may seem surprising that constraints away from the pivot wavenumber are finite given that 
a(n/ NI ) = oo, we remind the reader that the infinite running of /nl(&) is essentially multiplied with zero 
amplitude /^l when calculating the constraints at the fiducial value /nl(&) = 0. Closer inspection of 
Eqs. (4.3) and (4.4) confirms this argument. 



- 15 - 




10" 4 10" 3 10" 2 10" 1 10° 

k (h/Mpc) 

Figure 7: The same as Figure 6, but with a fiducial model /niX^) = 0. In the limit /j!j L — > there 
is no information on the running of non-Gaussianity ?i/ NL , and hence the LSS/BigBOSS constraints 
are sharply peaked and essentially constrain /nl(^) at only one wavenumbcr. 



we studied the extent to which a combination of LSS and CMB observations can constrain 
the running (Table 1) and /nl(&) as a whole (Figures 5, 6, and 7). 

For the power-law /nl(&)> we found that both the bispectrum measurement from the 
CMB Planck survey and power spectrum measurement from an LSS survey can constrain 
/nl(&) tightly in a relatively narrow range of wavenumbers around k ~ 0.1/iMpc . The 
scale best constrained by the CMB is larger (i.e. at a smaller k) than the scale best con- 
strained by LSS: we get complementary information about /nl(&) from the two data sets. 
The ability of LSS to constrain /nl(^) effectively at a wide range of scales depends on the 
survey parameters and the fiducial model of /nl(&) chosen, as is clear from Figures 5-7 
and Table 2. Nonetheless, large galaxy redshift surveys planned for the future may well 
be competitive with, or even better than, the constraints on the magnitude and running 
°f /nl(&) expected from Planck. 

Beyond the simple power-law model, we find that the combination of CMB and LSS 
helps pin down the best-constrained few principal components of /nl(^) better than either 
probe alone. Figure 4 shows that the degree of complementarity significantly depends on 
the details of (and systematics in) the LSS survey. 

The constraints from the DES and BigBOSS, and other upcoming LSS surveys can turn 
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out to be worse or Setter than those illustrated here, depending on how well the systematics 
can be controlled. While (for example) the photometric redshift errors [24], calibration 
errors [66], and assembly bias of galaxies [67] can all introduce parameter biases and degrade 
constraints, accurate calibration of these effects from simulations and observations, as well 
as selection of the "golden" class of objects with well understood properties whose clustering 
to use to measure non-Gaussianity, can cancel out these degradations. Moreover, we have 
not considered information from the LSS bispectrum which, while somewhat notoriously 
difficult to theoretically estimate due to non-Gaussian contributions from the gravitational 
collapse at late times (though see [68, 69] for recent progress on the matter), is nevertheless 
a very potent probe of primordial non-Gaussianity (e.g. [70, 71, 72, 73]). 

Overall, a full exploration of the LSS and CMB systematics is a herculean task beyond 
the scope of this paper; nevertheless, we think we captured a few key systematics with our 
choice of survey specifications and nuisance parameters. 

Finally, we introduced the figure of merit for measurements of non-Gaussianity, defined 
as the inverse area of the constraint region in the plane of non-Gaussian amplitude and 
running (see Eq. (4.5)). We are very encouraged by the fact that future constraints of 
non-Gaussianity will improve current-data figure of merit [64] by more than an order of 
magnitude, and thus shed interesting constraints on the physics of inflation. 
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A. Calculating the CMB bispectrum 

Calculating the CMB bispectrum is a problem that has been well-studied elsewhere in 
the literature, both for the general case and primordial local non-Gaussianity (e.g. [75]). 
Here, we briefly review the technique for calculating the bispectrum in the case of local 
non-Gaussianity, as well as the extension to the generalized local model that we discuss in 
this paper. 

The bispectrum is defined as: 

-^1^2^3,TJl"l2"l3 = Qf 2 m2 ^3 m 3 ) (-^-"l) 

where the ai m s are the coefficients on the spherical harmonic decomposition of the CMB 
sky. The a£ m s can be related to the Bardeen curvature perturbations $(k) by: 

a lm = J d 2 k ^y/ m (k) = ^{-if I $(k) 5 ,(fc)F4(k) (A.2) 
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Here, gg(k) is the CMB temperature radiation transfer function. There are several conven- 
tions used for this transfer function; ge(k) is related to the transfer function Ti{k) found 
in ([76]) by: 

gAk) = -. Tt(k) (A.3) 

Throughout this paper, we denote the radiation transfer functions as ti(k), defined as: 

t e (k) = -^gdk) = - 1 T e (k) (A.4) 

With these transfer functions, (A. 2) becomes: 

47r , f d 3 k 



^3 c&(k)^)Y/ m (k). (A.5) 



The angular-averaged bispectrum B^ 2 ^ 3 is related to the raw bispectrum Bi^ 2 i 3 ^ mim2rn3 
of (A.l) by the relation: 



B, 



hhta — 

mi,m2,m3 



2- \m im2 m 3 ) Be ^ mm - (A - 6) 



Here, (J^ ^ ^ ) is the Wigner 3j-symbol 6 . Substituting (A.l) and (A.5) into (A. 6), we 
obtain the following expression for the angular- averaged bispectrum: 



d 3 k x d 3 k 2 d 3 k 3 



rm m 2 m 3 J J (2n) 3 (2tt) 3 (2tt) 3 

xY/^koy/^^^ (a.7) 

Using the definition of the Bardeen curvature bispectrum, B&, 

($(ki)$(k 2 )*(k 3 )) = (2TT) 3 5(k 1 + k 2 + ks^fci, k 2 , k 3 ), (A.8) 

we find 

B ^i 3 = ^ E ™ 2 ? t 3 ) / ^ ^ <?**Y? imi (ki)y; 2m , 2 (k 2 )y; 3m , 3 (k 3 ) 



mi,m2,m:) 



xt <1 (fci)t <3 (fc 2 )t <3 (A;3)(y(ki + k2 + k3)B*(fci,fc2,fc 3 ). (A.9) 



(The prefactor of (—1)^1+^2+^3 vanished because the Wigner 3j-symbol ensures ^1+^2+^3 is 
even.) Taking advantage of several identities in [77] (their (12) and (13)), the orthogonality 
of the spherical harmonics, and the Gaunt integral identity ([5]), this becomes: 

Biihh = ( ~ ) hxMs I kfdki k^dk 2 kldk 3 B^(ki,k 2 ,k 3 )t ei (k 1 )t h (k 2 )t e3 (k 3 ) 



7T / 



00 

2 



x / r dr j ei (kir)ji 2 (k 2 r)ji 3 (k 3 r), (A.10) 







6 There are some computational difficulties that arise when evaluating the 3j-symbol for high £1,2,3; see 
Appendix D.2 for more on this. 
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where Ii x i 2 t 3 is the Gaunt integral 



^_ v 'M±a^±aek±a(^ft). (A . n , 

The real-space integral is now a one- dimensional integral in the spherical coordinate r, 
starting at our location and ending at infinity. This real-space coordinate is the difference 
in the conformal time An = f t ° ^ = c(tq — r e ) between the time when the CMB was 
emitted and the present. Nearly all of the contribution to the integral in r comes from 
a short period of time around the surface of last scattering, and there are no physical 
contributions beyond r > r max = tjq = ctq « 14.6 Gpc. To perform this integral, we 
sampled it 150 times between r max and r max — 2r*, where r max — r* is the comoving distance 
to the surface of last scattering. We also sampled 50 times between r max — 2r* and to 
capture any impact that late-time effects might have had. Increasing the sampling rate did 
not significantly improve our results. 

A.l Derivatives with respect to /nl and /nl(&) 

Using (A. 10) along with (3.3), we get the following expression for the angular- averaged 
CMB bispectrum in the constant /nl case: 

B ixiilz = 2A|/ NL (J^J j kjdh k\dk 2 k 2 dk 3 ^ ^^^^^^ + perm.^ 



oo 

2, 



x t ei (ki)ti 2 (k 2 )t£ 3 (k 3 ) / r dr ji 1 {k 1 r)ji 2 {k2r)j t3 {k- i r). 
Jo 

(A.12) 

Following [5, 53] we define functions ae(r) and /3e(r) to help us rewrite (A.12) as 



2 , ., 



7T 



a e (r) = - / k 2 t e (k)j e (kr)dk (A.13) 



7T 

Now Eq. (A.12) reads 



Mr) = 2 / k-( 2 - n °k e (k)j e (kr)dk. (A.14) 



/■oo 

B il£2 e 3 = 2Alf^I il £ 2 £ 3 r 2 dr (a ei (r)/3 e2 (r)Pt 3 (r) + perm.) (A.15) 

■/ o 

and hence 

o~r - -f— B hM 3 - (A.lb) 

C/nl /nl 

For the scale- dependent /nl(^) case, we use (3.5) to find that the angular- averaged 
CMB bispectrum is: 

= 2A}/ Waft / r 2 dr (^(^(^(r) + perm.) (A.17) 
01 /nl jo 



'NL 

where 



^uppcr 

aj(r) = - / 1 k 2 t e (k)j e (kr)dk. (A.18) 

7T y flower 
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A. 2 Polarization and cross-terms 



The bispectrum for multiple fields is a simple extension of the single field case. By analogy 
with Eqs. (A.l) and (A. 2), the multiple-field bispectrum is 

^W 2 <3,mim2m3 = ^Iimi a f 2 m2%m 3 )' (A. 19) 

where 

Air f d 3 k 

<n = y/W +T ) (-V' J (2^p ^)mYlM (A.20) 

and tUk) is either the temperature or polarization radiation transfer function. Using 
these definitions and running through Eqs. (A. 7) through (A. 17) again, we can rewrite the 
bispectrum for multiple fields if we just modify Eqs. (A. 13), (A. 14), and (A. 18) slightly: 

= \\ k 2 t p (k)j e (kr)dk; (A.21) 

0P( r ) = 1J k^ 2 - n °k p e (k)j e (kr)dk; (A.22) 
k u pp cr 



af{r) = - / k 2 t p (k)j e (kr)dk. (A.23) 
So for the constant /nl case, we have 

2 A 2 I^ lz ^ r 2 dr (a* (r)/3« (r)/% (r) + perm.) , (A.24) 



^3 



<9/nl 

while for the piecewise-constant /nl(&) case, we have: 

dB pqr r°° / n 

— ^ = 2 A^ 1 ^ 3 / r 2 dr (<(r)/3j 2 (r)^ 3 (r) + perm.) . (A.25) 

B. The covariance of the bispectrum 

It is usually a good assumption to consider only the Gaussian contribution to the covariance 
of the bispectrum, C. Using Wick's theorem, one can straightforwardly show ([78, 54, 55]): 

Ce 1 e 2 e 3 = C^C^C^ (B.l) 

where 

Cg = Cf V + <r|Wi = Cf v + Cf , (B.2) 

where C9 V is cosmic variance, while is the variance due to the noise and beam width 
in the survey; moreover, a 2 is the variance of the noise in the survey per pixel, and We is a 
"window" term relating to the survey beam type and width ([79, 80]). 7 For an experiment 



7 Note that [79] uses w 1 for what we are calling a 1 
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with multiple frequency channels (such as Planck or WMAP), the basic form of equation 
(B.2) still holds, but finding Cf is slightly trickier ([79]): 



1 



1 



£ 



1 



(B.3) 



V 



For uncorrelated Gaussian noise, o\{y) = cr 2 (v) is constant, and we can find its value for 
a particular experiment - for example, the Planck beam width and noise parameters are 
found in the Planck mission "blue book." 

We have only been dealing with temperature (TT), but it is not significantly harder 
to add in polarization (EE) and cross (TE) terms. The covariance matrix here is ([17, 54]) 



Noise is dealt with in the same way as in (B.2) for Cj T and Cf E in (B.5). Assuming 
that the noise for T and E are uncorrelated, a\ E = (ATAE) = (AT){AE) = 0, and thus 
Cf - TE = for all I. 

C. The high-peak limit 

Desjacques et al. [29] have identified a new term that contributes to the scale-dependent 
bias due to non-Gaussianity, which becomes important when the high-peak limit assump- 
tion is relaxed. This new term successfully explains previously mysterious discrepancies 
[10] between the theoretical expectation for the scale-dependent bias and the results of 
numerical simulations. Physically, the new term accounts for the scale-dependent mapping 
between the interval in the peak height dv (which is featured in the peak-background split 
derivation of the bias) and mass interval dM. 

Moreover, this term is only non-zero for cases when /nl ^ const, and therefore it 
affects constraints on /nl(&) that we study in this paper, but not the numerous forecasts 
for constant /nl studied previously in the literature. 

The new term corresponds to the second term of Eq. (2.3) 





(B.4) 



where 




(B.5) 



N(k) 



dlnF(k) 
d In ur 



(C.l) 



We can make the evaluation of this term more tractable by using the chain rule 




(C.2) 
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Now we will need to take the derivative of N(k) with respect to the f^ L , for our Fisher 
matrix. 



ON 



9fi 



NL 



F 



dan 
dM 



d 



dfhL 
d 



d<JR \ 

dM) df } 



i 

NL 



1 dF 

F(k)dM 

d f dF 



1 dF dF 
FdMdf^ 



(C.3) 



Equations (C.2) and (C.3) are everything we need to properly account for the new term 
in our Fisher matrix. Note that <tr and dan/dM are the only redshift-dependent quantities 
in N(k); since their redshift dependence is linear and exactly the same, it cancels entirely, 
leaving N(k) independent of z. 

The effect of this new term on the projected constraints for the f^ L , with a fiducial 
value of f^ L = 30, are seen in Figure 8. The figure illustrates that this new term removes 
much of the correlation between errors in neighboring and slightly broadens the range of 
scales at which the survey is sensitive to /nl(&)- Nevertheless, given that we are expanding 
our general /nl(&) model around a constant value (30 or zero), the effects of this new term 
on the constraints on the amplitude and running of /nl ~~ /nl an( ^ n /NL ~~ are sma H. 

D. Calculational Details 



D.l £ sampling and binning 

In evaluating equation (3.7), we do not actually use every 



< 



that would be incredibly 



computationally expensive. Instead, we sample and bin in I. We keep every £ up through 
£ = 40, at which point sampling drops off gradually until, at t > 100, only every tenth £ is 
sampled. The "width" of the bins in £ are given by the equation 

1 [{ii - £_x) + (£ i+1 - £i)} = - £ 1 . 1 ). (D.l) 



A, 



D.2 Calculating the Wigner 3j-symbol 

We need to be able to calculate the Wigner 3j-symbol for large (> 1000) values of £1,2,3 hi 
order to evaluate many of the expressions we're interested in. Unfortunately, the 3j function 
built in to the GNU Scientific Library can't properly evaluate the symbol for £1$ 3 > 70. 
Thus, we were forced to create our own special-purpose 3j-evaluator. Thankfully, we're only 
interested in the special case 771^2,3 = 0; as it turns out, in this case, the 3j-symbol reduces 
to (see Wolfram Mathworld: http://mathworld.wolfram.com/Wigner3j-Symbol.html): 



h £ 2 £ 3 




-if 



(2g-2£ 1 y.(2g-2£ 2 )\(2g-2£ 3 y. 



(2g + iy 



{g - £^.(3 - £ 2 )\{g - £ 3 )\ 



if L = 2g; 







if L = 2g + 1, 
(D.2) 

where L = £\ + £ 2 + £3. Since (D.2) involves evaluating the factorials of relatively large 
numbers when any of £1,2,3 are large, we used Stirling's approximation to perform the 
factorials - but we needed the factorials to remain accurate even when the arguments were 
small, so we used six terms in the approximation. 
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fc(h/Mpc) Mh/Mpc) 
(c) Unmarginalized with the new term (d) Marginalized with the new term 

Figure 8: Illustration of how the inclusion of the correction to the scale-dependent bias from 
[29] affects the forecasted constraints on the /j5j L from DES. For comparison, the green line is the 
constraint found for a constant /nl using the same assumptions. 
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